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Abstract 



We show that taking into account the interplay between perturbative and non-perturbative 
effects, the power-suppressed shift to the broadening distributions becomes B dependent, and 
the non-perturbative contribution to the mean values becomes proportional to 1/ {Q\Ja s {Q)). 
The new theoretical treatment greatly improves the consistency of the phenomenology with 
the notion of the universality of confinement effects in jet shapes. 
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1 Introduction 



There is mounting evidence in favour of the universal pattern of leading 1/Q non-perturbative 
power corrections to collinear /infrared safe (CIS) observables 0. They include jet-shape 
observables such as the thrust T, the C-parameter, squared jet masses and the jet broadening 
B (definitions can be found, for example, in 0). 1/Q power effects are expected also in the 
energy-energy correlation function and in e + e~ annihilation. The 1/Q power terms are 
seen both in the means and the distributions. For instance, the distribution of events in T 
in the 2-jet region 1 3> 1 — T ^> A/Q was predicted [||,|4]] to be given by a simple shift of 
the corresponding perturbative QCD prediction by a constant value inversely proportional 
to Q. The same pattern was expected to hold for other jet shapes V as well, with the true 
distribution in the kinematical region 1 ^> V ^> A/Q being related to its perturbative (pt) 
counterpart by a shift, 
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with cy an observable dependent, perturbatively calculable, numerical coefficient which for 
the thrust, C-parameter and total and heavy-jet masses, e.g., is 
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The parameter V oc 1/Q is the non-perturbative (np) quantity which effectively "measures" 
the intensity of the QCD interaction over the infrared momentum region. Its magnitude 
can be interpreted as being related to a mean value of the QCD coupling over the infrared 
region, say, k < fij = 2 GeV, 



dk a s {k) . 



1.3) 



A proper definition of V includes an infrared matching scale \xi which is necessary for merg- 
ing, in a renormalon-free manner, the pt and np contributions to a given observable. It 
also calls for a two-loop analysis of the non-perturbative contribution || , without which the 
magnitude of the power correction cannot be quantified better than up to a factor of order 
unity. 

Experimentally, np effects in the thrust distribution have been found to be consistent 
with the shift rule ( |1 . 1| ) , with cto — 0.5. The same value was experimentally extracted 
from the Q-dependence of (1 — T) 0-0. The C-parameter has also been studied (both the 
distribution and the mean) and found to be consistent, with a similar value of «o Pf- 

Power effects in the broadening B were also expected to shift the distribution, but with 
Ab logarithmically enhanced (A# cx In Q/Q) ||. The HI collaboration however stated that 
the data was not consistent with a InQ enhancement 0. Most recently the JADE collab- 
oration studied the discrepancy and showed that relative to the perturbative distribution, 
the experimental distribution is not only shifted, but also "squeezed" ||10|| . 

In this paper we revisit the broadening. We show that the coefficient of the power cor- 
rection shift, QAb, is neither proportional to InQ, nor a constant, but rather is a function 
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logarithmically depending on B. The reason is the following. The non-perturbative contri- 
bution to V comes from the emission of gluers (gluons with finite transverse momenta, of 



the order of the QCD scale A, with respect to the quark direction ||11|| ). For instance, in 
the soft approximation, the contributions to the thrust, C-parameter and broadening from 
a secondary parton i can be expressed as 

(l-n = ^e-^ } Q = -^-, 2B l = ^, (1.4) 

ictf cosn T\i \v) 

where the transverse momentum ku and the rapidity T]i are measured with respect to the 
thrust axis. The feature that 1 — T and C have in common is that the dominant non- 
perturbative contribution to these and similar shapes is determined by the radiation of soft 
gluers at large angles. This radiation is insensitive to a tiny mismatch, (Q q ) = 0(a s ), 
between the quark and thrust axis directions which is due to perturbative gluon radiation. 
Therefore the quark momentum direction can be identified with the thrust axis in the pt 
and np analysis of T and C. 

The broadening, on the contrary, accumulates contributions which do not depend on 
rapidity, so that the mismatch between the quark and the thrust axis could be important. 
As we shall see later, for the pt analysis of broadening this mismatch is irrelevant, to next- 
to-leading accuracy. However, the mismatch plays a crucial role for the np effects in the 
broadening, both in the means and the distributions. 

If one naively assumes that the quark direction can be approximated by that of the thrust 
axis, as is the case for the perturbative contribution, then B accumulates np contributions 
from gluons with rapidities up to the kinematically allowed value rji < f] maiX ~ \n(Q/k t i). In 
this case one would find the shift in the B spectrum to be logarithmically enhanced, 

A B = Cb P (1.5) 

where cb = 1(|) for the total (single-jet; wide-jet) broadening 0. What this overlooks is 
the fact that the uniform distribution in rji (defined with respect to the thrust axis) holds 
only for gluon rapidities not exceeding | ln0 9 |. Hard gluons with energies k^i > ku/Q q are 
collinear to the quark direction rather than to that of the thrust axis and therefore do not 
contribute essentially to B. As a result, for a given quark angle, the np contribution SB to the 
broadening B comes out proportional to the quark rapidity (we note the distinction between 
the np contribution 5B to a given perturbative configuration, and the np shift Ab to the 
perturbative distribution, where the latter is integrated over all perturbative configurations 
leading to particular value of B). For the single-jet broadening one has 

SB ~ cxV-ln—. (1.6) 



Let us describe, semi-quantitatively, how ( |1.6|) affects the np corrections to (B) and to the 
B distribution. The power correction to the mean single-jet broadening (B) 1 is obtained by 
evaluating the perturbative average of 5Bi in (|1.6| ), 



<B>P = <5> 1 -<*>r ) ^^.(lni->. (1.7) 



2 



The PT-distribution in Q q at the Born level is singular at Q q = 0. In high orders this 
singularity is damped by the double-logarithmic Sudakov form factor. As a result, the np 
component of (B) 1 gets enhanced by 

1 \ 7T / \ 

In— ) ~ — . 1.8 

For the mean wide-jet broadening (B) w the result has the same structure with the replace- 
ment Cp — > 2Cp due to the fact that now it is radiation off two jets which determines the 
Qq distribution. 

The shift in the single-jet (or wide-jet) broadening can be expressed as 



A 1 (B) ~ Cl V-(]n±-J , (1.9) 



where the average is taken over the perturbative distribution in the quark angle Q q while 
keeping the value of B fixed. Since Q q is kinematically proportional to B, the log-enhancement 
of the shift in the B spectrum becomes 

A 1 (B) ~ Cl V-\n^-, (1.10) 

with B = B Q (a s \nB) a calculable function slowly depending on B. Thus, the shift in the 
£>i (B w ) distribution becomes logarithmically dependent on B. 

Finally, the shift in the total two-jet broadening distribution At{B) will also be derived. 
It has a somewhat more complicated B dependence. In the kinematical region where the 
multiplicity of gluon radiation is small, a s In 2 B <C 1, one of the two jets is responsible for the 
whole PT-component of the event broadening, while the second is "empty". That "empty" 
jet contributes the most to the shift: in the absence of perturbative radiation the direction 
of the quark momentum in this jet stays closer to the thrust axis. This results in 

A T (B) ~ A a (2J) + <B>i NP) ~ c{P (in I + ^== + O (1)) . (1.11) 

In these circumstances the B dependence of the total shift practically coincides with that of 
a single jet. In the opposite regime of well developed pt radiation, a s In 2 B>1, the jets are 
forced to share B equally, and we have 

A T (B) ~ 2- Ai(B/2) ~ 2-dV ln-^- (1-12) 

B 



In the present paper we derive the correspondingly modified predictions for the B de- 
pendent 1/Q shifts in the broadening spectra and the 1/Q corrections to the means, which 
supersede the earlier results of @. The paper is organised as follows. 

Section 2 introduces the quantities to be studied, and recollects the next-to-leading- 
logarithmic perturbative result for the single-jet, wide-jet and total distributions. 

In section 3 we derive the new results for the power corrections to the means and distri- 
butions. 
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Section 4 presents a comparison of these new results with experimental data. 

In the concluding section || we discuss our results in the context of the interplay between 
pt and np effects. 

The remainder of the paper consists of appendices which contain technical details of 
the derivation and a Monte Carlo study which illustrates some important features of the 
analytical results. In the last Appendix [F] we give the full list of formulas to be used to 
describe the leading power correction to the broadening means and distributions. 



2 Broadening distribution 



2.1 Distributions and means 

Broadening in e + e~ annihilation with c.m.s. energy Q is defined as a sum of transverse 
momenta of final particles with respect to the thrust axis, 

2Q-B X = £|£u| . (2.1) 

The symbol X here marks the set of selected final particles. The three options under dis- 
cussion are 

1. single-jet broadening, Bi, with the sum running over particles in one hemisphere (right, 
B 1 = B R , or left, B X = B L ), 

2. total broadening of the event Bt = Br + Bl, and 

3. wide-jet broadening which is the larger of the two, Bw = max{BR, Bl} on an event 
by event basis. 

We define the integrated single-jet distribution as 

™ = S(B) S § = (2.2) 

For the purposes of the present paper it suffices to treat the double differential distribution 
in Br, Bl as the product of two independent jets, 

i d 2 a i da -, da . . „ . . . „ . . 

o" — — = o" — a' 1 - — = ^[(Br) ■ ^[(Bl) . (2.3) 

dBRdBL dBR dBL 

In this approximation the total and the wide broadening distributions are related to Si 
through 



dB 
dZ w (B) 



dB 



I dB 1 E' 1 (B 1 )E' 1 (B-B l ), (2.4) 
Jo 

2-£ / 1 ( J B) / dBN^BN); E W (B) = E?(S) . (2.5) 
Jo 
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Single-jet distribution. The Mellin image of the single-jet B distribution, a{y), is defined 



as 

where the contour runs parallel to the imaginary axis (in the second case we place the contour 
at Re v > in order to satisfy the normalisation condition (|2.8|) , see below). u{v) is a regular 
function of v in the entire complex plane. It is limited (decreases) at Re^ — > +00 which 
results in 

E[{B < 0) = £1(5 < 0) =0. 
In the left half-plane it increases exponentially, 

a{v) oc e - 2Bm \ Rezv^ -00, (2.7) 

with B m the maximal value of single-jet broadening. This ensures the kinematical con- 
straints, 

MB > B rn ) = a(0) = 1 , (2.8) 

E[(B > B m ) = 0. (2.9) 



Total broadening distribution. The Mellin representation for the integrated total broad- 
ening distribution in the factorisation approximation follows from 



E T (B) = / a» . (2.10) 



Mean single-jet broadening. To calculate (B) we average B with the differential distri- 
bution to write 

(B) 1 = / BdBE[(B) . (2.11) 
Jo 

The property (|2.9| ) allows us to extend the B integration in ( [2.1 1| ) to infinity, provided that 
the ^-contour has been shifted to run to the left of the imaginary axis, Re^ < 0: 

f°° , r c+l °° dv 2R „ , x r c+ico dv aUA . „ N 



Mean total broadening. The representation (|2.10|) immediately leads to 

d 



= -o-'(o) =2(B) 1 . (2.13) 

u=0 
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Mean wide-jet broadening. The mean B\y, according to (|2.5|) , is given by the integral 
BdB ^f- = -J dB {2\Sh(B) - 1] + - l] 2 } 

Jl +o / ^77777^ (2-14) 



2m v 2 2 J 2m v\ [y + v\ 

Rcz/<0 I, Re^i<0 



f dv g(u)[l + a(-u)} [™ d" <u)(t(-u) 

J 2m 2v 2 ' 2^W + 




2m 2 v 2 



The latter integral cannot be evaluated by residues, since the integrand exponentially in- 
creases in both directions, o{y)o{—v) oc exp(| Reu\B m ). 



2.2 Kinematics of small B 

Here we introduce the kinematical variables needed to analyse B in the soft limit. We use 
the Sudakov representation for the momenta of the primary quark and antiquark p\ and P2 
and of the accompanying soft partons hi 

p 1 =A 1 P + B 1 P + Pa , 

p 2 = A 2 P + B 2 P + pa , (2.15) 
fci = aiP + (3iP + k ti . 

The two light-like Sudakov vectors P and P are taken along the thrust axis and their sum 
is the total incoming momentum (2PP = Q 2 ). 

In the soft region all quantities 

an, Pi, (I- A,), (1-B 2 ),^ 

are small and of the same order while A 2 and B\ are much smaller (quadratic in transverse 
momentum) and will be neglected. As shown in Appendix 0, for the two-loop perturbative 
contribution one can approximate the quark- antiquark direction with that of the thrust 
axis, i.e. one can neglect pa and pt 2 . However, as explained above, the non-perturbative 
contribution is sensitive to the difference between the thrust and quark axis and then we 
shall take into account pu- We introduce the quark and antiquark angular variables (rescaled 
quark transverse momenta) 

P = -r-ptl, p = —p t2 . (2.16) 

A x t> 2 
The soft radiation matrix element depends on the invariant 

7 \ — - - 7T ( k u ~ ocip) 2 {ku - Pip} 2 , 

where the first (second) factor is the squared transverse momentum of fcj with respect to the 
quark (antiquark) direction. 
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In the following we concentrate on soft partons emitted within the right hemisphere, 
near the quark direction. The kinematical constraint for such a parton is oti > fy. Since a 
massless emitted parton satisfies a^i = k 2 JQ 2 , 

a i >k ti /Q. (2.17) 

For such a right-hemisphere parton one can neglect the difference between the thrust and 
antiquark directions (hi — fop ~ k ti ), so that the matrix element depends only on the 
transverse momentum with respect to the quark axis 

km = hi - a iP- ( 2 - 18 ) 

According to the definition of the thrust axis, the vector sum of the parton transverse 
momenta in the right (or left) hemisphere is zero and we can write 

= pti + Yl = a ip+ -p+Y Ru ■ ( 2 - 19 ) 

ieR ieR ieR 

where we have neglected the small a-components of the left-jet partons (1 — A 1 ~ J2en a i)- 
The right-jet broadening B R is defined by 

2B R Q = p a + hi = A lP + J2 k ti ^ P + - a iP ) , (2.20) 

ieR i&R ieR, 

where again we have neglected the small a-components of the left-jet partons. 
Similarly one introduces Bl, the left broadening, 

2B L Q ~ p + J^iku-fop). (2.21) 

i<=L 



2.3 Resummation and the radiator 

Here we recall the resummed expression for the broadening distribution at small B [|EJ. In 
this region the distribution is obtained by resumming contributions from the emission of any 
number of soft partons and can be expressed in terms of a "radiator" . The two-loop analysis 
of the radiator is presented in detail in Appendix [A]. 

The integrated single-jet broadening distribution (e.g. the right jet) is given by 

S 1 (5) = Q 2 I ^-5 2 (p + Y^)® [2BQ-p-J2(ku-a lP ) \ , (2.22) 

n ieR V ieR J 



where we have used ( |2.19|) to represent the transverse momentum constraint in terms of the 



rescaled quark momentum p ( |2.16| ) and the secondary parton transverse momenta with respect 
to the quark direction, K ti . The factors da n /a are the n-parton emission distributions which 
factorise for B < 1, i.e. in the soft limit. To perform the sum over n one needs to factorise 
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also the constraints. This is obtained by using the Fourier and Mellin representations for 
the delta and theta functions 



6 2 (p + *h) el2BQ-p- - cup) 

i&R \ i&R / 

= f ^ b e ~ibp f dl> e 2Bu e ~up/Q "Q e -v( ktl -a tP )/Q e -ibK U _ (2.23) 

J (2tt) J 2ixiv i( _ R 

where we have introduced the vector impact parameter b conjugate to the transverse mo- 
mentum with respect to the quark axis ( |2.18| ). One can now resum soft-parton emissions in 
the Mellin representation ( [2.6|) . The resulting integrated single-jet broadening distribution 
is 

Si (B) = [ ^-e^ B a(u) , *(v) = [ e -W e -» P /Q e -R(*,M (2 . 2 4) 

J 2m v J (2tc) 2 

The exponent R in ( [2.24j ) is the "radiator" analysed in Appendix |A]. Since soft real and 
virtual partons are responsible also for the non-perturbative effects, the radiator contains, 
together with the pt contribution, also a np correction, 

R(i/, b; p) = + TZ^ . (2.25) 



PT component. At two-loop accuracy, the pt part of the radiator is p- independent. At 



large values of v and/or bQ, it is a function of the combination |L2 



?1E 



v + Jv 2 + (bQ) 2 



and reads 



= n{p) = jf* f n'(Q/k) , n'(p) = ^a s (Q/p) (W - ^ . (2.27) 



(2.26) 



The running coupling is taken in the physical scheme |H| . This expression accommodates all 
terms a™P. n+1 and a™£ n with £ = In p. We present 1Z as a function of a single dimensionless 
parameter p but it is implied that it also contains a Q-dependence via the running coupling 
a B (k t ), with Q/p < k t < Q. 



NP component. According to the procedure developed in 14 , the leading power-suppressed 



np contribution to the radiator has the following general form, see Appendix |A.2j for details 



k .10 



dm 2 



Sa e ff(m 2 ) ■ fl (m 2 



(2.28) 



Here M. is the Milan factor which emerges from the two-loop analysis, and 8a e fi(m 2 ) is the 
np component of the effective coupling related to the standard a s by the dispersive relation, 
with m 2 the corresponding dispersive variable acting as the gluer's "mass" . 
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The factor Q (m 2 ) is a "trigger function" which is specific to a given observable. A 
power-behaved np contribution is determined by the leading non-analytic in m 2 — > term in 
Qq. For the broadening measure, as well as for other jet shapes, the leading non-analytic 
piece, 50o, is proportional to \Ar? and is given by the following expression, 

5n (m 2 ) = is-— ■ / — / — (Jl + u 2 + 2ucos^ -u) . (2.29) 

Q 7p/q u J 2n \ ) 

Here 

w = / = — e ' , U < 77 < In — 



Q m 
with the gluer rapidity with respect to the thrust axis. The expression in the brackets in 
( p.29| ) accounts for the mismatch between the quark and thrust axes. Since, with account of 
the azimuthal integration, the integrand falls rapidly at u — > 00, the rapidity-integral, in the 
m = limit, has a finite value proportional to In p. Thus the np component of the radiator 
becomes 

ft( NP ) = v-Vhi— , r/ = In ^ - -0.6137056. (2.30) 

v Q 



The definition of the NP-parameter V oc l/Q introduced in || is recalled in Appendices [A. 2 
and 0. 

We notice that if in (2.29) we formally set p = 0, which corresponds to disregarding the 
perturbative quark recoil, we reconstruct the old wrong answer p| 

KW^U * ^ • - |) , (2.31) 

where we have restored the 3/4 "hard correction" coming from the region a ~ 1. 

In the following we recall the perturbative part of the B distribution, while the non- 
perturbative one will be treated in the next section. 



2.4 Recollection of the perturbative result 



The p-integration giving the perturbative part of the single-jet distribution cr(u) in Q2.24j) is 
explicitly carried out in Appendix [A.l|, see (|A.28|), 



d 2 pd 2 b 



1 v 



H = e lB v 



1+y 



(2.32) 



To evaluate the Mellin integral we introduce the following convenient operator representation 
which we shall use extensively below, 



a=0 



(2.33) 



The Mellin image of the integrated single-jet distribution can then be represented as 



(y) = e -K{e- d «) ( e ^\(a))- a ■ v- a 



a=0 



(2.34) 
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where we have introduced the function |T2 

• i 



A" 



i y 



dy fl + y 

2 



dz 



o 



Performing the //-integration we arrive at 

dv 



Using the identity 



2iixv 



7 2Bv ^(PT) 



V) 



1 + z 
2z 



1 



A(0) = 2, A(oo) 
'e 7E A(a) 



1 . 



x g{a) 

we can absorb the power factor x~ a , 

X(7Z')e^ 



a=0 



r 1 + a 



2B 



(2.35) 



(2.36) 



a=0 



a=0 



X 



n' = n'i 



dK{x) 



2B r/ln.r ' 

into rescaling of the argument of the perturbative radiator: 



-Tl(xe- 9a ) 



X(a) 



(2.37) 



(2.38) 



a=0 



T(l + a) 

Performing the logarithmic expansion of the radiator, 

-71 (xe~ da ) = -7l(x) + 7l'(x)d a - \7l"{x)d 2 a + . . . , 

we conclude that the action of the operator on a function which is regular in the origin 
reduces to substituting 71' (x) for a, while 71" (x) = O (a s ) and higher derivatives produce 
negligible corrections: 

e- n{x) g(7Z'(x)) (1 + O (7Z"(x))) . (2.39) 



a=0 



Thus we can evaluate ( |2.36|) as 



Ef T) (5) 



;i+o(« s )) 



(2.40) 



T(l+7l'(x)) 

where the parameter x is a function of B implicitly defined by ( |2.37| ). This expression can 
be simplified. To this end we observe that within our accuracy we can substitute B^ 1 for x 
in 71' 



71' (x) =7l'(B- 1 ) + 0(a^ +1 \n n B) , n>0; K '(B' 1 ) 



2C F 



7T 



a s {QB) InB- 1 



At the same time, the finite product xB should be kept in the exponent 7Z(x). We finally 
obtain 

Ef T) (S) = — , B = , ft' = ft'(£T 1 ). (2.41) 



r(l + 71') ' s ' e^X{TZ') ' 
For the wide-jet and total broadening distributions we obtain, in a similar way, 

e -2TC(B- 1 ) e -2n{B- v ) 



X W (B) 



MB) 



(2.42) 



r 2 (l + ft') ' ±y y T(l + 27l') 

These results are equivalent to those based on the steepest descent evaluation presented 



in 



12]j . The resummed pt answers can be given in various forms which are equivalent at the 
level of the leading a™ln n+1 B and next-to-leading a™ ln n B contributions to the exponent, 
see, e.g. equations (4.25) and (4.26) of ||I2"| . The specific form which should be used for 
matching with the exact second order pt answer will be discussed in Appendix [F]. 
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3 Non-perturbative correction to the B distribution 



With account of the leading np contribution the integrated single-jet B distribution takes 
the form (cf. ( EQjp ) 



jT bdbe~ n ^ ] jf V d P e- vp/Q (^j J (bp). (3.1) 
Performing the p-integration, see ( A.39| ) in Appendix |A.2 , the distribution assumes the form 

tT ( I /) = crCPT)^) + + ; (3 2) 

with the non-perturbative correction given by the following expression 



We recall that 



dy 



IE ~ VO 



In 



i + y 

2z/y 2 




(3.3) 



l p 7i5 
2 e 



(bQY 



v e 



IE 



1 + y 



(3.4) 



The function f(u) will determine the np corrections to the B distributions and to the means. 
In particular, the np effects in the distributions are given by the inverse Mellin transform of 
f(u) which we calculate in terms of the operator technique introduced in ( |2.33| ). In the case 
of the pt distribution, the operator exp{— 7Z(xe~ 9a )} acted upon a function g(a) regular at 
a = 0, see ( |2.39| ). In the case of B T , as we shall see later, we need to calculate the operator 
exp{— 7Z(e~ 9a )} acting upon a function g(a) oc 1/a which is singular in the origin. To this 



end we introduce and evaluate in Appendix |B.2| the function 



E(x) 



D TZ(x) n(xe- a <>) -1 



a=0 



dz 



(3.5) 



The value E(l) enters into the expression for the power correction to (B), 



E(i) = = i 3 ^ 

1 ; 2^Cf~cT s 4 6C F 



+ O 



'OL* 



(3.6) 



As far as the pt component of the distribution er( PT )(z/) in ( |3.2j ) is concerned, we have derived 
the resummed expression which is valid for large values of v and determines the B distribu- 
tions in the B 1 two-jet kinematics (soft limit). The B distributions at large B, as well 
as the means (B) , are determined by finite moments v. The soft resummation programme is 
irrelevant here, and the exact fixed-order analyses should be carried out instead. At the same 
time, the np component j[y) in ( |3.2j ) which originates from soft gluers remains applicable 
both for B <^ 1 distributions and the means, that is for arbitrary values of v. Having said 
that we turn to the calculation of the leading np effects in mean broadenings. 
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3.1 Power corrections to means 



Substituting the sum of perturbative and non-perturbative components, ( |3.2| ), into the ex- 
pressions for the mean broadening ([2.12|) and ( |2.13|) we obtain 

(B) T - (B)^ = 2 ((S) a - (B)™) = -V /(0) . (3.7) 

The value of / in the origin is related with the function E ( |3.5| ) in Appendix |B], 

/(0) = -(E(l) + 77o) [1 + 0(a s )\. (3.8) 

This results in 

^-^ PT, ^(vfe + i-^ + "°)' (3 - 9) 



Mean wide-jet broadening is given in terms of the Mellin integral in (|2.14| ). Extracting 
the non-perturbative component we obtain 

(B} w = (B) ( l T) +V[5-±f(0)} +0{V 2 ), (3.10) 

with 

6 = \ fZ 2^ VM*™^ -/HV PT V)] • (3.H) 



This integral is evaluated in Appendix resulting in 

S = l(^ e - 2 ^( e - ett )-e-^ e " Sa ))a- 1 | (3.12) 

We note that the first operator here corresponds to E(l) evaluated with 2TZ substituted for 
the radiator. This leads to 



V 

2 

V( v 3_ ft \ 

2 V2^U 4 12C F 7,0 J 



(B) w - (B)^ T) = V - (e-W-^a- 1 + Vo) [ 1 + O (a s ) ] 

" " (3.13) 



The result can be obtained from that for the single-jet (total) broadening by the simple 
substitution CV — > 2Cf in (B.9I). 



3.2 NP shift in the single- and wide-jet distributions 

Having obtained the 1/Q power corrections to the means we now turn to the analysis of the 
np effects in the distributions in the B « 1 region. To this end we need to perform the 
inverse Mellin transform of the distribution ( |3.2j ) containing the np term f{y) given in ( |3.3| ). 
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To perform the integration over y (related to the impact parameter b, see ( |3.3|) ) we make 
use of the operator representation (|2.33|) to write 



dy 



i y 2 



1 + y 



1 + y 

2y 2 



2 - Je - Vo + hi + d a -y) ■ v~ 



a=0 



The result of the y-integration can be represented as follows: 

f{v) = e~ n ^ (e^A)" a (^2 - lE - In 2 - m + p{a) - X {a) + d a -^j v~ a 
where A = A (a) is defined in ( |2.35|) and p and x are given by the integrals 



a=0 



p(a) 



XW 



dy (1 + y 



1 2/ 2 V 2A 



In 



— + — 



dy ( 1 + y 



y V 2A 



- = 2 v-i 

a 



(3.14) 



The functions p(a) and x{ a ) are regular at a = and vary slowly between p(0) = — 2 + 21n2, 
p(oc) = In 2 and x(0) = 2 In 2, x(oo) = 2 correspondingly Introducing the function 

C(a) = 2 - 7b - In 2 - r? + p(a) - x(a) , 
the operator expression for the np component of the distribution takes the form 



e^\)- a \C{a) + d a --)u- a 



(3.15) 



a=0 



The ^-integration of Q3.2|) is now readily performed. For the single-jet distribution we have 



E 1 (B) = e- nie ' da \e' E \)- a 



^- [ C{a) + d a -~ 



a J Via) 



(2B) a 



T(l + a) 2B 

We observe that the potential singularity at a = cancels in the combination 

(2B) a 



(3.16) 



a=0 



a / r(a) 



T(l + a) 



D(B,a) 



with 



—D(B, a) = C(a)+ln2 J B-<0(l + a) = lnS - 770 + 2 + p(a) - x(a) + + a) 

a function regular in the origin, a = 0. The function ip(x) is defined as i/j(x) = d\nT(x)/dx. 
We finally arrive at 



E^B) = e 



-■R(e- 9 ") 



e-*E\(g) 
2B 



T(l + a) 



(W)Y 

r(i + a) 



aV ~ 
l-^D(B,a) 



l-^D(B,a) 



a=0 



a=0 



(3.17) 
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Here we have absorbed the power factor x~ a with x defined in ( [2. 37j ) into a rescaling of 
the argument of the pt radiator as we did before. As in the case of the pt distribution 
considered above, the action of the operator on a regular function results in substituting 
Tl'(x) ~ TZ'iB- 1 ) for a, according to ( ggg) . We get 



E 1 (B) 



e -K(B~ l ) 



(3.18) 



with B defined in ( 2.41]) . This correction can be cast as a -B-dependent 1/Q shift of the 
perturbative distribution, namely 



E 1 (B) ~ Y>r ] [ B - -D.iB) 



V 



(3.19) 



with 



d^b) = d(b,k') = \iiB- 1 + r] -2- p(n') + x (n') + ^(i + n')-^(i), n' = n'(B- 1 ). 

(3.20) 

For not too small values of B such that 1Z' -C 1 the shift -Di assumes a simple form 

£>! ~ lnS^+jyo, ft'<l, (3.21) 
while in the opposite limit of extremely small B, 



Dl ~ ln 2B + lE + m ' U ' > 1 



(3.22) 



The same shift applies also to the wide-]et broadening distribution which, to the needed 
accuracy, is simply given by the squared single-jet distribution according to (|2.5|). 



3.3 NP shift in the total broadening distribution 



To obtain the integrated distribution for the total broadening, T,t(B), we need to perform 
the inverse Mellin transform of o 2 [y) in (|2.10j) . Invoking the operator representations ( |2.34|) 
and ( 3.15 ) for and f[y) correspondingly, we construct the product (a, 6^0) 



a 2 M 



-n{e-^)-n{e- d b) (lE 



(e lE \(a)Y a {e' E \{b)Y 



1 + 2-uVi C{a) + d a -- 



-(a+b) 



+ O (V 2 ) 



a=b=0 



where we have made use of the a <-> b symmetry. Evaluating the i/-integral we obtain 
E T (B) -- < 



{2B)a+b I 2 V (C(a) I fl M 
r(l + a + 6) +2 2B{ C{a) + da ajT(a + b)\ 



a=b=0 
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Here the first term gives £( PT ) (5) and the second one accounts for the leading 1/Q correction 
contribution. Recalling the definition of x fl2.37Q we write (7Z' = TZ'(x)) 



M B ) 



e -K(xe- a * ) e -1Z{xe- 9 b ) \ A ( a ) 



T(l + a + b) 

l + ^T(l + a + 6) [C(a)+\n2B + d a -- . 

B \ a J 1 (a + b) 



(3.23) 



a=b=0 



The first factor is identical to that for the perturbative distribution. The second factor can 
be represented as 



B 



(a + b) (C(a) + In 2B - + a + b)) - - 

a 

V 



l-_((o + 6) D(B,a)+ip(l + a + b)-ijj(l 



- . (3.24) 



We observe that, in contrast to the single-jet case (cf. ( |3.17| )), we get a correction term 
singular in a, which reminds us of the case of (B) considered above. With the unity and the 
non-singular term in ( |3.24|) we proceed as before substituting VJ for a and b, see (|2.39|) . The 
factor a + b produces 2TZ' which makes it possible to interpret the correction in terms of a 
shift. The contribution of the singular piece in ( |3.24j ) to the shift is calculated in Appendix |D[ 
The final result reads 



MB) 



,(PT) 



B - ^D T (B) 



D T (B) = 2D X (B) + 2[V>(1 + 2K') - ^(1 + K')) + H(B^ 1 ) , K' = TZ{B~ l ) 



(3.25a) 
(3.25b) 



Here the single-jet shift Di(B) is given in ( |3.20 ), and B defined in ( |2.41|) . The np shift in 
the total broadening distribution, ( ^.25| ), includes the function HiJB^ 1 ) which is analysed in 
Appendix [D[ The shift has a rather complicated B dependence: it changes from 



D T {B) 



to 



In — + const 
B 



1 



for ^WB<1 



D T (B) w 2 ■ ln^ , for— In 2 > 1 

B ix 

Indeed, for moderately small values of B such that ct s In 2 B < 1 (1Z' < 
(see (03)) 

ix 3 p 



/a^), we have 



H(B~ 1 ) 



hxB 



2VC^~ S 4 6C F 



Together with (|3.21|) this results in 

D T {B) ~ 2 ■ (ln^- 1 + 770) + H^B- 1 ) ~ D 1 (B) + (D) . (3.26) 
Here the first term is the single-jet shift, D\ ~ ln(po/(BQ)) = \riB~ 1 + t]q; the second term, 

A) 



(D) 



7X 



3 



2^CpCi, 4 6C F 



(3.27) 
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is nothing but the np correction to the mean single-jet broadening (|3.9| ). 



" (B)^ = j (D) ■ (3-28) 

The physical origin of this result can be simply understood. In the region under consid- 
eration, (TZ' < y/oQ, the multiplicity of perturbative gluon radiation is small. The non- 
perturbative shift is then dominated by the fluctuations in which one of the two jets is re- 
sponsible for the perturbative component of the event broadening: Bl <C Br « B^, or vice 
versa. Non-perturbative effects in the narrower jet are better pronounced: in the absence 
of perturbative radiation the direction of the quark momentum stays closer to the thrust 
axis thus bringing in a large contribution (D) oc \/^/oT s to the shift D T . It describes the 
non-perturbative correction to the mean broadening of the narrower jet. This contribution is 
practically B independent: its residual B dependence emerges at the level of a O (^ln 2 B) 
correction. In these circumstances the B dependence of the total shift coincides with that 
of a single jet. 

In the opposite regime, a s In 2 B > 1, multiple perturbative radiation becomes a way of 
life. Here both jets take responsibility for perturbative broadening, and the coefficient of 
InB in D? grows bigger than unity. It doubles in the limit of extremely small B values, 
a s ln B>1 (TZ' ^> ^ya^), where the differential spectrum flattens off (and then starts to 
decrease, when TZ' > 1) because of severe Sudakov suppression. In this regime the jets are 
forced to share B equally, and we should expect 

D T {B) ~2£>i(£/2)- [1 + ^((asln 2 ^)- 1 )] , TZ' » y/a£ . (3.29) 

Indeed, at a s \n 2 B > 1 (IZ' > ^/a^) we have H(x) ~ l/TZ' <C hi-B -1 . Near the maximum 
of the distribution, TZ' > 1, the relative size of the if -contribution becomes as small as 
H/D T ~ (TZ')' 1 /InB- 1 = O (a B ). The first two terms in (|3.25b|) then combine into fl3~29l). 



4 Comparison with experimental data 



4.1 Means 



Using formulas flF.12| , |FTT3"D we perform fits to the data for the mean total and wide-jet 
broadenings and compare to the values obtained for the thrust, the heavy-jet mass and the 
C-parameter (with fixed-order perturbative coefficients taken from M). 



The results that we obtain are (fitting to data from |6|, |7|, |l5l- 19| and references in [20 
partially based on a pre-existing compilation |2T|j) 



Variable 


a s 


a 


X 2 /d.o.f. 


Bj< 


0.1170 ±0.0023 


0.4508 ±0.0225 


14.9/(23 - 2) 


B w 


0.1189 ±0.0025 


0.3911 ±0.0305 


12.8/(22 - 2) 


1 — T 


0.1177 ±0.0013 


0.4976 ±0.0087 


57.0/(40 - 2) 


C 


0.1206 ±0.0021 


0.4527 ±0.0110 


10.7/(10-2) 


Ml/Q 2 


0.1171 ±0.0012 


0.5602 ±0.0224 


15.2/(27-2) 
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• DELPHI 

■ ALEPH 

T L3 

X JADE 




Figure 1: Energy dependence of (B) T and (B) w . The dotted line is the LO result; the 
dashed line is LO+NLO and the solid line, LO+NLO+power correction, (|3.7| ) and ( |3.13| ) 
respectively. 



Figure [I] shows the a comparison between the data and the fits for B? and B^. 
4.2 Distributions 

We consider only the Bt distribution, for the sake of illustration. The relevant equations 
are ( |F.4| , [F.5| ) for the non-perturbative shift and the perturbative spectrum is as discussed in 



section |F.2| . We fit 8 data sets (using log-i? matching, with the distribution constrained to 
go to zero at the kinematic limit) [f?|,p!6 HT9i|2^| and obtain 





a 


X 2 /d.o.f. 


0.1158 ±0.0007 


0.5368 ±0.0077 


68.7/(59 - 2) 



Two examples of the distributions are shown in figure 0, together with JADE (35 GeV) and 
OPAL (91.2 GeV) data. 



5 Conclusions/Discussion 

In this paper we have demonstrated that the sensitivity of the broadening measure to large 
parton rapidities has made the non-perturbative shift in the B spectra B dependent. As a 
result, the pt distribution gets shifted to larger B values and squeezed on the way, in accord 
with the recent experimental findings || . This result supersedes the previous expectation of 



17 



the In Q-enhanced shift [[| . It is worthwhile noticing that the correct result contains a single 
universal NP-parameter V which puts broadening on an equal footing with other jet shapes 
1— T, C, etc. The previous (incorrect) result contained an additional NP-parameter, Qb in 
(|1.5|), related to the log- moment of the coupling defined similarly to ( |1.3| ) but with an extra 
factor In k under the integral. 
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Figure 2: B T distributions compared to JADE (35 GeV) [f£J and OPAL (91.2 GeV) |T§ 
data. The dashed lines are resummed distributions without the power correction, while the 
solid lines are the resummed predictions with the power correction. 
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The reason why the shift in the B spectra has become B dependent is the interplay 
between pt and np gluon radiation effects. The radiation of gluons with fixed transverse 
momentum does not depend on the gluon rapidity as defined with respect to the quark 
direction. At the same time, the broadening accumulates the moduli of transverse momenta 
of final partons with respect to the thrust axis. The latter, starts to differ from the quark 
direction when the normal pt radiation is taken into account. This mismatch does not 
matter, at next-to-leading accuracy, for the pt analysis but plays a crucial role for the np 
effects both in the B means and distributions. 

If one naively assumes that the quark direction coincides with that of the thrust axis then 
the contributions from gluons (gluers [|TI]]) with finite transverse momenta, k t ~ Aqcd ; and 
rapidities up to the kinematically allowed value 77, < r] max ~ \n(Q/k ti ) sum up to provide the 
log Q-enhanced np shift in the B spectrum, see ( [2.31| ). However, gluers with large energies, 



koi ^ k ti /Q q are collinear to the quark direction rather than to that of the thrust axis and 
therefore do not contribute to B because of the CIS nature of the observable (gluons collinear 
to the quark increase B by exactly the same amount by which the quark contribution to 
B is reduced due to the longitudinal momentum recoil). As a result, the np contribution 
to B comes out proportional to the quark rapidity, logl/G g , with Q q the quark angle with 
respect to the thrust axis which is due to radiation of pt gluons. Since, kinematically, 
Qq oc B, averaging logl/O g over the pt distribution in g results in the log-B enhancement 
of the np shift. 



NP effects in the presence of PT radiation. Disregarding the ever-present pt radiation 
is known to produce confusing results in a number of cases. For example, the first-order 
(one-gluon) analysis of the np effects in the heavy-jet squared mass produced the wrong 
expectation: adding a single gluer to the qq system in e + e~ (as the third and only secondary 
parton) and constructing M 2 of the qg system we find a 1/Q confinement contribution to the 
squared mass of the heavy jet, the one our gluer belongs to. Meanwhile, the opposite lighter 
jet acquires neither a pt nor a np contribution to the mass. As a result the np correction to 
Mfj comes out equal to that for thrust, 

ct = c M 2 = c M 2 r , c M 2 = . (5.1) 

In reality there are always normal pt gluons in the game which are responsible for the bulk 
of the jet mass: Mjj/Q 2 ~ a s > M 2 /Q 2 ~ a 2 3> SM^ P . In these circumstances it is the pt 
radiation to determine which of the two jets is heavier. The gluer(s) contribute equally to 
both jets, 

ct = c M 2 = 2c M 2 h = 2c M 2 . (5.2) 

It is worthwhile remarking that experimental analyses carried out before 1998 were based 
on the wrong expectation ( |5.1|) . 

Another important example of the interplay between np and pt effects is given by higher 
moments of jet shapes, e.g. ((1— T) n ). For such a quantity one obtains, symbolically, 

(( (1 _ r)pT+ ^J), a , +Q ,^ + ...(^J, „>, 
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The leading power-suppressed contribution is at the level of 1/Q. At the same time, the 
one-gluer analysis for such an observable would formally produce 



only 



VII 



■n 



which np correction is suppressed as a high power of 1/Q. In the presence of pt radiation, 
however, the leading 1/Q contribution is still here, though reduced by the a s (Q) factor but 
far more important than the 1/Q n term. 

Another "mistake" of this sort brings us closer to the B issue. Consider the transverse 
momentum broadening of the current-fragmentation jet in Deep Inelastic Scattering (DIS), 
that is the sum of moduli of transverse momenta of particles in the current jet. Adding a 
gluer to the Born (parton model) quark scattering picture we get three equal contributions 
to B: two contributions from the quark p which recoils against the gluer k emitted either in 
the initial (IS) or in the final (FS) state, \p±\ = \k±\, and one contribution from the gluer 
itself when it belongs to FS. Taking into account the pt radiation, however, the FS quark has 
already got a non-zero transverse momentum, ~ a s ■ Q, a substantial amount compared 
to kj_ ~ Aqcd- In this environment the direct gluer's contribution is the only one to survive: 
the np recoil upon the quark gets degraded down to a 1/Q 2 effect after the azimuthal average 
is performed, 



The true magnitude of the 1/Q contribution turns out to be a factor three smaller than 
that extracted from the one-gluer analysis. 

A preliminary health report. Now that the loophole in the theoretical treatment of 
the jet broadenings has been eliminated, one can return to the much debated issue of the 
universality of confinement effects in event shapes. 

To illustrate the overall consistency of the universality hypothesis we show in Figure ||] 
2-standard-deviation contours in the a s -c*o plane for a range of means and distributions. 
The curves labelled "old" are the results to the fits || using the old (wrong) formulas. 
The situation with the total broadening distribution is greatly improved by the updated 
theoretical treatment. We expect the wide broadening distribution to be equally improved, 
but this remains to be verified. 

The fits for a s and «o from the mean values are also generally consistent with each other 
and with those from the distributions. However, the agreement between different event 
shapes is still not perfect. In the case of the heavy-jet mass we believe that this may in part 
be related to the treatment of particle masses, which have more effect on jet masses than on 
the thrust or the C-parameter (which are both defined exclusively in terms of 3-momenta). 
We leave this question for future consideration. 

Outlook. Another important issue is that of the power correction to the jet-broadening 
in DIS. Formally the extension of our results to the DIS case is quite a non-trivial exercise. 
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Figure 3: 95% CL contours for jet shape means (dashed) and some distributions (solid). The 
curves for the T, C and old Bt and Bw distributions are taken from || . The curves for the 
means are to be taken as purely indicative since we have not accounted for the correlations 
between systematic errors (which, where available, are added in quadrature to the statistical 
errors). Additionally for some observables we may not have found all the available data. 

As a first step it would be necessary to carry out a resummed pt calculation for the DIS 
broadening. This has so far not been done. 

The situation for the mean broadening measured with respect to the thrust axis is fairly 
simple though, since (modulo factors of two associated with the definition of the broadening 
in DIS f2~3fl ) it is equivalent to a single hemisphere in e + e~: 

<£>Dis, thrus t - = v ( - n + \-^ + ^ + (yffi) ■ ^ 

\2^C F acMw{Q) 4 bLj F J 

For the mean broadening defined with respect to the photon (z) axis the situation is more 
complicated because of the dependence on perturbative initial-state radiation. To a first 
approximation, at moderate x, one can view the DIS event as a rotated e + e~ event where 
the broadening is measured in the right hemisphere with respect to the axis of the quark 
in the left hemisphere: i.e. the relevant transverse momentum for determining the rapidity 
available to the np correction is p = \pi + P2\- Since this is very similar to max(pi,p2) we 
have a situation like that for the wide-jet broadening, and the leading power correction is 
suppressed by factor \[2 compared to (|5.3|) : 

(b w - <*>£E - * { ^mkm + J " ife + * + ° ' (5 ' 4) 

Even though we have chosen to include some sub leading terms of O (1), it is likely that there 
are other terms of 0(1), arising through the x dependence of the problem. In particular, 
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for the important case of small x which is presently being studies at HERA the process is 
dominated by the boson-gluon fusion mechanism and the analogy with the two quark jets in 
e+e- gets lost. Parton multiplication in the initial state of the DIS process is more intensive 
(t-channel gluon exchange) and increases with In 1/x leading to increasing characteristic 
transverse momentum of the struck quark, which contributes to -Bdis,2- 

Still more damaging will be x-dependent effects at 1 — x <C 1 where phase-space restriction 
on PTtransverse parton momenta may destroy the leading term through a dependence 

on ln(l — x). 

So (|5.4j ) can be used only for moderate x, and even then one should remember that it is 
subject to subleading corrections of O (1). 

Finally let us remark that 1/Q power effects can be envisaged (and should be studied) 
also for other jet characteristics such as shape variables of 3-jet events (thrust minor, acopla- 
narity) as well as, for example, in the distributions of the accompanying Et flow in hadronic 
collisions and DIS. 
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A Perturbative and non-perturbative radiator 



The two-loop radiator in the soft limit is given by (see 

da d 2 K t ( a s (0) 



R(u,b;p)=AC F 



a hk 



An 



+ X {kI) [1 - u(*)] 



+ 4CV J dT 2 {h,k 2 ) (^) 2 ±M 2 (k 1} k2) [l-u(h)u(k 2 )}, (A.l) 
where, for the right-broadening distribution, the source u(k) is given by (see ( [2.23 )) 



u(k) 



1 _ e -v(kt-ap) IQ e ~ibR t 



${a - kt/Q) 



k t = K t + ap 



(A.2) 



with the right-hemisphere constraint ( j2.17|) included. 



The function x{ K t ) i s ^ e virtual correction to one-gluon emission. In the physical scheme 
which defines the coupling as the intensity of soft-gluon radiation [13|, it can be written in 
terms of the dispersive integral 



dfi 2 k\ 



o ^ 2 {i4 + ^ 2 ) v 4vr 



-2C A ln 



(A.3) 
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The collinear divergence in x is compensated by that of the real two-parton contribution 
described by the matrix element M(ki,k 2 ). Taken together with an ill-defined a s (0) of 
one-gluon emission, they participate in forming a finite running coupling. 

The two-parton phase space in ( |A.1| ) is 

dT 2 (k u k 2 ) = dm 2 — ■ dz -- , (A.4) 

71 a 2,71 

The first three variables a = ct\ + a 2 , Kt = Kf\ + K t2 and m 2 = (ki + k 2 ) 2 are those of the 
parent gluon k, while z, 1—z and <fr are the momentum fractions and the relative azimuth of 
the two secondary partons, qq or gluons: 

— * — * 

cti = za , a 2 = (1 — z)a , m 2 = z(l — z)q 2 , q t = . (A. 5) 

Hereafter we choose <fi as the angle between 2- vectors qt and Kf 

The probing functions u(ki), u(k 2 ) depend on all the parton variables. In order to extract 
the contribution responsible for the running coupling we split 

l-uihMh) = [ 1 — u(k) ] + [ u(k) — u(kx)u(k 2 ) ] , (A.6) 

where we have introduced a "probing function" u(k) for the parent gluon. There are various 
ways to define the source u(k) if k is massive. The prescription we choose consists in 
substituting k, 2 — > n 2 + m 2 in the massless expression (|A.2|) . In particular, the transverse 
momentum with respect to the thrust axis, k t , 

k 2 t = k 2 + (ap) 2 + 2apK t cosipkp , 

gets replaced by 

k' t 2 = k 2 + m 2 + (ap) 2 + 2ap\J k 2 + m 2 cos^fep • (A. 7) 

Thus we define the inclusive probe (first term in (|A.6|) ) by 

1 - u{k) = (l - e~ v{k 't- ap)IQ e- ih V* 2 t+ m2 cos *hb\ . #( a _ k'JQ) . (A.8) 

The inclusive contribution depends only on the parent gluon variables, so that the two-parton 
matrix element M 2 (k\, k 2 ) can be integrated over z and <fi to give (see 0) 

dz ^M 2 (k u k 2 ) = 1 {-A, + 2C A ln ^M±Hfl\ , (A.9) 
27r2! m^im 1 + nf) { m 4 J 

with p = ^Ca — \nf the first beta function coefficient. We treat the (3q term of this 
equation together with a s (0) to form the so-called naive contribution, the logarithmic term 
together with x( K ?) {inclusive correction), and the second term of the trigger function ( |A.6| ) 
as the non-inclusive correction. 

At the perturbative level, the naive term provides the dominant contribution, while the 
inclusive and non-inclusive terms stay at the level of the next-to-next-to-leading two-loop 
correction. As far as the non-perturbative 1/Q correction to the perturbative radiator is 
concerned, all three give comparable contributions, the latter two providing the so-called 
rescaling Milan factor to the naive one. 
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Naive contribution. The naive contribution reads 



Ro(f,M = 4C F 



\°^8{m*) - A ■ n (n 2 + m 2 ) . (A.10) 



Here we have introduced the "naive trigger function" representing the [1 — u(k)} factor 
integrated over a and ifj^p- 



Q (k? + m 2 ) 



da 



2tt 



1 - e 



-u{k' t -ap)/Q -ibyj K^+m? cos Vfcf, 



(A.ll) 



with fcj defined in (|A.7 ). The lower limit of the logarithmic a-integration, corresponding to 
separation between the two hemispheres, is actually k' t /Q. However for p <C Q, this limit 
can be approximated by the p-independent value given here. 



Inclusive correction. The inclusive correction can be represented in terms of a difference 
between the naive trigger functions for m^fl (real) and m = (virtual contribution) as 

, , x ^ ^ f dm 2 /a s \ 2 f dn 2 , K 2 (K 2 + m 2 ) . , 9 9n ^ . 9nN 

iuMjp) = sc^cu J—r{-£) J ^rf-i in A * m4 (fioK 2 + m 2 ) - fi (* 2 )) . 

(A.12) 



Non-inclusive correction. Finally, the non-inclusive correction describes the mismatch 
between the actual contribution to Br from two partons and that of their parent: 

R ni (z/, b;p) = 4C F J dm 2 dn 2 ^ dz ^M 2 (m 2 , k 2 , z, 0) ■ Q ni , (A.13) 

with the non-inclusive "trigger function" 

Q ni = C d _^l d ± [u{kl + h) _ U (A; 1 ) W (A; 2 ) ] . (A.14) 
Jo a J 2n 



A.l PT contribution to the radiator 



In the perturbative evaluation of the naive contribution we neglect m 2 C K( C Q 2 in the 
trigger function and use the dispersive relation for a s , 

r^A{ <mm2) -^} = <A - i5) 

to obtain 

/■"» r J 2 

Ri PT) (^6;p) = — / • fioK 2 ) • (A.16) 

y k 2 

The inclusive and the non-inclusive terms do not contribute to the pt radiator at two loops. 
This is due to the fact that, in spite of the singular behaviour of the matrix element, M 2 oc 
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1/m 2 , the m 2 -integrals in ( |A.12| ) and ( A.13|) converge, because the trigger functions Qo(k 2 + 
m 2 ) — Q (k 2 ) and VL ni vanish in the collinear parton limit. As a result, these contributions are 
determined by the non- logarithmic integration regions K t ^m^Qjv and provide negligible 
corrections to the radiator of the order of a 2 ■ In v a 2 In B , with a single-logarithmic 
enhancement factor originating from the a- (rapidity-) integration. 

The pt expression (|A.16|) can be greatly simplified, to our accuracy, by observing that 
one can neglect the difference between quark and thrust axis and set p = in the trigger 
function (|A.11| ). Moreover, for large v and b one can substitute the trigger function Qq by a 



suitably chosen cutoff in the phase space integration. 



Negligible p ^ effect. As we shall see later, the effects of quark recoil (p 7^ 0) modifies 
the structure of the answer for the non-perturbative contribution. However, at the perturba- 
tive] level a small departure of the quark direction from that of the thrust axis proves to be 
negligible, producing a O (a 2 In B) subleading correction which lies beyond the scope of the 
approximation, i.e. keeping under control terms O (a"ln m B) with m >n [|I"2"1 ,|2^ 



To verify this we observe that the difference between the pt radiator (|A.16|) , (|A.11|) and 
its p = value, 

R( PT >(M;0) = 9L [ Q d -^a s {K t ) ln£ [ 1 - J (b Kt ) ] , (A.17) 

7T Jo Kf K t L 

in the large-^ limit is a function of the ratios vp/Q and bQ/u, 

L( PT )(6,,;p)-R( PT )(6,,;0) = / (^,^) + O 0) . (A.18) 



To verify this we introduce the rescaled dimensionless variables q = ^ and u = bQ/v to 
write 

PT PT C F f 1 da _ f v d 2 q f 1 - e ~<*(<i~»p/Q) e -i^(q-up/Q) 



2 



2vr Jo a Jo tt { {q-vplQ) 2 q 

Eq. ( |A.18| ) follows from the fact that the g-integral converges and can be extended to run 
up to q = oo instead of q — v. As a result, it does dot depend explicitly on v except 
through the combinations vpjQ and bQ/v. In the essential region, bQ ~ v ~ Q/p ^> 1, 
these combinations are of the order one, so that 

/ vp bQ 



Radiator as a function of a single variable. 

v and b is obtained by the replacement ^1 



[\-e-^/ Q Jo{bk t )] =>Q(kt-Q/ji) 



The final simplification of ( |A.17| ) for large 
+ V^TTbQj 2 



/i = e 



IE 



(A.19) 
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which introduces a negligible error O (a s ) to the radiator and gives, in the next-to-leading 
logarithmic approximation, 



(A.20) 



Here we have replaced the soft da /a factor in the trigger function (|A.11| ) by the exact q — > qg 
splitting probability to account for a "hard" subleading correction originating from the region 
a ~ 1. The result is the —3/4 term in the integrand. While necessary for achieving single- 
logarithmic accuracy in the perturbative treatment, the hard part of the kernel is irrelevant, 
as we shall see below, for the non-perturbative contribution which is dominated by the region 
a ~ K t ~ m -C Q. 

We give expressions for the radiator for three forms of the coupling: 

• For fixed coupling 

%£),£ = hip. (A.21) 



7l 



Taking into account the one-loop running of a s we obtain 



^ (PT) (/i) 



4CV 



o Jo 



d£' 



4C F 



L — - 

4 



lnll- r 



+ 



To evaluate the radiator with the two-loop running coupling we use 



to obtain 



Here 



a 8 (Q) 



T = 



2tt 



0o(L + j\nL) 



7 



0i 



Pi 



19 

51 - -n f 



In?. (A.22) 
A v ' 



(A.23) 



AC F 
~0o~ 



2tt 



(T-l) ln--(T-T, 

-L i 



7 



T- 



bin 2 



T 

Tn. 



2tt 



0o ot s {Qln) 



+ 7 



(A.24) 



(A.25) 



0o a s (Q) 

and a s is the running coupling in the so-called physical scheme |13| which can be related to 
the standard Ms coupling by ( |b'.7| ). 

The radiators with the two-loop and one-loop a s deviate at the level of a O (af In 3 
term which contribution is under control and should be kept in the pt distributions. At 
the same time, for the sake of simplicity we use the (two-loop) radiator with the one-loop 
coupling ( |A.22| ) for evaluating the np contributions to the means and spectra. For example, 
in the means we know that O appears as a constant term, suppressed with respect to the 
leading term by a factor of yfa^, see ( |3.9| ). It is clear therefore that a 0\ contribution will 
enter at the level of a term of O (a s ), which is beyond our control. 

In conclusion, the pt part of the Mellin transform of the single-jet broadening in ( |2.24| ) 
reads 



[ bdbe- n(FT) ^ [ pdp e~ up/Q J (bp) . (A.26) 
Jo Jo 
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The p-integration gives 







Q 2 u 2 y 3 v 



which results in 



a^\v)= r^e-^W; ^=\e^{l + y)v. (A.28) 

Ji y 



A. 2 Non-perturbative contribution of the radiator 

To extract the power-suppressed contribution to the radiator we use the procedure developed 
in Ul4| . It is based on introducing the effective coupling a eff (m 2 ) related to the usual coupling 



constant a s via the dispersive relation 



Qs{k) r dm \ a fjl . (A.29) 



k 2 J (m 2 + k 2 ) 2 

We then substitute the non-perturbative "effective coupling modification" 5a e s for a e s and 
look for the leading non-analytic in m 2 term in the m 2 — > limit. 

To obtain the np contribution to the radiator we substitute 
in the naive, inclusive and non-inclusive contributions. 



Naive contribution. We obtain 

n„ r Q2 ( -H \ r Ql Jk 2 



-.2 , „ 



— / ^ a eS (m 2 ) n (m 2 ) + O (a s (Q)) 



(A.31) 



At the pt level, this expression is equivalent to (|A.16 ), with the non- logarithmic perturbative 



a s (Q) correction coming from the region m 2 ~ Q 2 



To trigger the leading power correction in ( |A.31|) w e substitute 5a e s for a e s and consider 



the leading non-analytic in m 2 term VL{m 2 ) oc yrn? which is obtained by linearising Vt in 
m ~ Kt, Qa <C Q. In this approximation the np component of the trigger function, SQ, does 
not depend on b, and we get 



2 rvl^fl du ,2k ^ 



5n (m 2 )=u-—- — / — V1+u 2 + 2ucos^-m , (A.32) 

Q Jp/Q u J 2n \ ) 

where we have introduced u = ocp/ y/rn? as an integration variable. The ^-integral converges 
and is determined by the region p/Q<u< 1. Therefore, in the m 2 — > limit we can replace 
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the actual upper limit, p/ym? ^> 1, by oo (neglecting the 0{m 2 /p) contribution to <5il, 
which is analytic in m? and thus does not produce a np correction). We have then 



Q 



pip), p = v[\n^ + 0[^ 



Po 



P 



P 



(A.33) 



where the integration constant po/Q is given by 

2^ 



In 5 



Q 



00 du 
u 







(•\/ 1 + u 2 4~ 2m cos ^ — u) — $(1 — it) 



-0.6137056 = 77o 



Recalling the definition of the first non-analytic moment of 5a e s, 

• m 5a eS (m 2 ) , 



2tt 



dm 2 



o 



we finally obtain 



TV 



dm 2 2 2 2Ai po 

— — oa e fj(m ) ■ os 2(m J = v ■ In ■ 







p 



(A.34) 



Milan factor. It is straightforward to verify that the inclusive and non-inclusive trigger 
functions, in the linear approximation in n t ~ m, are proportional to the same function p(jp) 
that determines the naive trigger function oT2 given in ( |A.33| ). We have 

SQ in = p- ^k 2 + m 2 - kA , (A.35) 
SVl ni = p ■ (kh + K t2 - a/ k| + J . (A.36) 

Such a structure is typical for l/Q power corrections to various jet shapes and leads to the 
universal rescaling of the naive contribution ( |A.34| ) by the so-called Milan factor, for details 



sec 



With account of the Milan factor, the full two-loop np component of the broadening 
radiator reads 

7^ Np ) = Tl^-M = v-Vln^, V = ^L. (A.37) 

p Q 



In conclusion, the Mellin transform of the single-jet broadening in (|3.1| ) reads 

(r{y) = j bdbe~ n{fl) J pdpe- up/Q J Q {bp)(^j . (A.38) 
We need to evaluate the p, b and v integrals. The result of the p integration reads 



pdp e -»p/Q il j o{b p) = -l^lf- , y = v (A . 39) 

o \PoJ v y v 

where I<p(y,y) is related to Legendre function and has the following small- V expansion, 

lp(».u) = V(2 + ,'V)[uy%) yPi + w(-) 

Q J \V' (A.40) 

l + vv(2- lE -r) + ln^-\nv-y) + O (V 2 ) ; r, = ln^. 
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B Calculation of /(0) 



To calculate the non-perturbative correction to (B) 1 = \ (B) T we need to find the value of 
/ in the origin. To this end we split f(y) given in (3.3) into two pieces: 

f ,s f°° vhdh (n ,i v+VK r Tl?\ _ nip) 

h{u) = l WT^w[ 2 - lE - m + ln ^WT^V) e ' (B.i) 

M) Jo V + v* 6 ■ 
Here and in the rest of this section, for simplicity, we measure b in units of Q. 

After extracting — In v from the first term, the remaining integral in y = \Jb 2 + v 2 / v 
converges and is therefore determined by the region b ~ v — > where the radiator can be 
dropped as O (a s ). We get 

fi(y^0) = 2- lE -r] + p(0)-\n2u + O (a s ) ~ In 2 - V o - 1e ~ In v . (B.2) 

The second piece which is logarithmic in b we integrate by parts to write 



f 2 (u - 0) = -\ In (b 2 + v 2 ) e- n W\ h b Zo + \ [°° d (e~ n ^) In (b 2 + v 2 ) 

Jo 

POO 

= lnu+ d (e- n(fi) ) In 6 , p, -> ±e 7B 6 , with v -»• . 



(B.3) 



In the region of finite 6 = (9 (1) the true radiator is (a s ) and its contribution to the integral 
there can be disregarded. To make it more explicit we cut the integral from below at some 
value b > bo ~ 1 and integrate by parts back again: 

db 



poo PC 

/ d{e- n ^ Eb ^)lnb = -\nbo-(l + 0(a s ))- / 

Jbo Jbo 



e -Tl{e'Ebl2) 



The answer to the required accuracy does not depend on bo, which can therefore be chosen 
arbitrarily. It is convenient to take bo = 2e _TB in order to represent 

f°° d7 

/ 2 (i/-0) = ]nu-Qn2- 7E ) - / — e "*M , (B.4) 

which finally gives 

/"°° 

/(0) = -Vo- - e- n{z) = -Vo ~ E(l) . (B.5) 



B.I Mean straight from the distribution (consistency check) 



We used the residue of a{y)jv 2 at v = to calculate (B) according to fl2.12| ). Here we verify 



that the same answer follows, within our accuracy, from integrating the spectrum: 

dB B — (E X (B) - 1) = - J dB (E X (S) - 1) 

dB (p l (B)-E? ,T \B)). 



(B.6) 
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The np component of the mean can be obtained directly by integrating (|3.16 ) over B up to 
the kinematical boundary B m ~ 1: 

,-1 /9R \ a 

(B.7) 



(B) - (B) (PT) = V c -K(e- a «) C ( a ) ~ ^i 1 + Q ) + H 2B r, 



r(i + a) 

Setting a ~ 7?.' ~ a s to zero in the non-singular pieces we get 



(PT) 



V 

' 2 



ln^-^o-e-^) 



■ a] l ( B m) a [Maj 



(B.8) 



a T(l + a)e a ^ 

In the final term, singular in a, the last factor should be expanded to first order in a. Using 

(a5o) = 1 + ° ^ ' r(1 + a)ga7B = 1 + ° ^ ' 

we observe that the dependence on the kinematical boundary value _B m cancels, 



and the result given in (|3.9|) is reproduced. 



ln^ + e-^)- 



B.2 Calculation of E(x) 



We define the function 



E(x) 



e n(x) e -n(xe- 9 ") a -i 



a=0 



°° dy a 



(B.9) 



(B.10) 



Replacing 

1 

a Ji y 

we apply the rule ( |2.39| ) to absorb the y~ a factor into a rescaling of the argument of 7Z, get 
rid of the differential operator and represent E in terms of the logarithmic integral, 



E(x) 



dy 

i y 



-Tl{xy) 



dz 



Jl{x)-K{z) 



(B.11) 



The perturbative radiator (as a function of a single variable x > 1) is defined in ( |2.27| ). 
The first coefficients of the logarithmic expansion of the PT-radiator with the one-loop a s 
(£ = lnx, L = InQ/A) are 



W {x) 
K"(x) 
TZ'"(x) 

K {n) (x) 



(in a; - |) 



2C F a s (Q/x) a s (Q/x) 



AC F t-\ 2C F a s (Q/x) 

Po L-l n 
4C F L-\ 

A> (L-l) 2 ~ 

8C F L — | f3p 

& (L-£) 3 " 2CV 1 ; ' a s (Q/x) ' 

fn-D! //5 2 a s (Q)' 



7T 



a s (Q) 



»\2 



(2TZ"(x)V 



n-2 
2 



16vrC F 



(B.12) 
(B.13) 
(B.14) 

(B.15) 
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with Q = Qe 4 . We also point out the structure of the combination 

TZ'(x) 



C F a s (Q) 



71 



(lux - f) 



The next step is to substitute (lnx — d a ) for lnx to construct the operator 
K{x) - TZ(xe- 9a ) = TZ'(x)d a - \n"{x)d 2 a + \n'"{x)dl + . . 
We start by noticing that 



F(a;K") = exp {-\K" d 2 a ) a' 1 



71 



N 



2-Jl" Vv 7 ^ 77 



where the function N is related with the probability integral, 



N(t) 



2e 



tr r°° 



dxe~ x = e* (1 - $(*)) = -=e l Erfc(t) 

V 7T 



and has the following behaviour: 



2t o 4£ 3 



JV(f) 



1 



7lf 



1 3 

+ ^ + 



£ < 1, 
t > 1 



As a result, 

F(a; ft") ~ a -1 for a > \/W' . 
Now we introduce the first derivative to obtain 

exp {K'd a - \Tl" d 2 a ) a' 1 = e n ' d " F{a; 11") = F{a + K'; K") 

fa + K r 



71 



N 



2TZ" V V2W 7 

To estimate contributions of higher derivatives, n > 3, we use ( |B.15| ) to derive 

^ d n N(t) 



n(n) {da)n -F{a + K'; 11") - ^ " s(0) ^° QsiQ) 



nl 



lQCpn a s (Q/x) \ 16ttCf 



dt T 



t=s+- 



(B.16) 



;B.17) 



;b.is) 



(B.19) 

(B.20) 
fB.21) 



(B.22) 



(B.23) 



We conclude that n = 3 contributes at the level of 0(1), while the contributions of higher 
derivatives are down by the factor {^/a^) n ~ 3 <C 1. 



Evaluating the third derivative, 

F(a + 1Z , ;1Z") 



exp 



K"'dl 



6 



1 + ^^)1 (a + 7?': 7?"' 



we obtain 



1Z"'d 3 

^-^F(a + 1Z';1Z"] 
o 



K 



a=0 



3 (ft 



a s (Q/x] 



X(s) 



[B.24) 



(B.25) 
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where the function X has been introduced, 

^d 3 N(t) 



X(t) = 

It has the asymptotic behaviour 

X(t) = 
X(t) = 
For ( |B.9|) we finally obtain 

E{x) = 



8 dt 3 



zdz e - z ~ 2t ^~ z . 



(B.26) 



l-^t + 4t 2 + 0{t 3 ) , t«l 



4t 4 



+ o (r 6 ) , t > 1 . 



(B.27) 
(B.28) 
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;b.29) 



The neglected terms in (|B.29| ) are of relative order a s . 



For s = TZ'/V^R 77 < 1 (see flETT^ )) we substitute for iV and X the expansions ( |B~20| ), 
( p.27|) and obtain, keeping contributions up to O (a/c^), 



E{x) 



7T 



2s 



2K" 

2^/U F ~a B {Q/x) 



+ s< + ... 



7T 



6C, 



(lnx 



+ 



\fC F a s 



(lnx 



3n2 _ Po 

" ' A> 6C F AC F 



^/C F a s (\nx - I) 



2a/CV 



7T 



a s (Q/x) +/ 2 (lnX ^ 



/3 y/C F oi s 



+ 



(lnx 



3\2 



We get 



7T 



3 

mx H — 



0o 



+ O (v^ln 2 x) . 



For x = 1 (the operator entering the expressions for mean broadening) this gives 



(B.30) 



:B.3r 



E(l) 
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00 



2^/C F a s (Q) 4 6C F 



;b.32) 



C Calculation of £ 



Here we calculate the Mellin integral necessary to determine (B) w , 

1 r°° dv 



2 7-ioo 2vriz/ 



[/0V PT ^)-/(-*V PT V)] 



(C.l) 
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Invoking the operator representations Q2.34 ) and (|3.15|) for and f(v) correspond- 

ingly, we obtain 



S = e -K(ze-da) e -K{ze- a b) 



with z = 2e lE and 



dv / / \ / t \ \ / Ma) 



X(b) 



F{a) = C{a) +d a --. 

a 

It is implied that we have to set a = b = after applying the differentiations. Now we 
introduce the real integration variable v = —iv, add the negative and positive f-beams 
into 21m, and start the f-integration from some finite value vq = 0(1) so as to ensure 
applicability of the \aige-v logarithmic expression for the pt radiator. This gives 

-6 



S = e -K(ze- d <>) e -Tl(ze- a b) 



1 °° — (F(a) - F(b)) v-^ sin < h Z*l(M)~ (M 

VQ 2ttv 2 V 2 / V 2 



Observing that the regular pieces C(a), C(b) cancel in the difference F(a) —F(b) at the level 
of O (a s ), and that the ratios A/2 produce negligible corrections O (a 2 + b 2 ), we are left with 



a + b 



2tt 



a + o I a b 



a=b=0 

+ 0{{a-bf) 



a=b=0 



where we have absorbed the power of vq into additional rescaling of the arguments of the pt 
radiators 1Z. Evaluating the derivatives and using the a <-> b symmetry we arrive at 

1 1' 



8 



1 -1l(v ze- a *) -H{v ze- a b) 



a + b a 



a=b=0 



1 ( -2-n(v \e-^) _ -n(v Xe-^)\ -1 



a=0 



The finite rescaling of the argument of the 1Z operator by the factor vqz produces, according 
to ( |B.31| ), a subleading correction \d.(vqz) = O (1) which is of the relative order ^fal and 
should be kept under control. These subleading corrections however are identical for the two 
terms and cancel in the difference, thus ensuring independence of the result on the arbitrary 
parameter v , at the O (a s ) level. We conclude, 



8 



a=0 



(C.2) 



According to ( |3.10| ), to obtain {B) w we have to add to (C.2) half of the np correction to 
single-jet (B), that is — 1/(0) with /(0) given in ( p.5| ). The main piece of the latter cancels 
the subtraction contribution described by the single-jet operator applied to 1/a, and we 
finally arrive at 



(B) w - <£>£ T) 



1 V ( e -2-K{e-«° ) a -l 



7T 



a=0 

3 

+ - 



00 



(C.3) 



2^2C F a s {Q) 4 12C F 



This result can be obtained from that for a single-jet (total) broadening by the simple 
substitution Cf — > 2Cf in ( |B.31|) . 
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D Analysis of 



Consider the singular piece of ( |3.24j ), 



S = e -n(xe- d -) e -Tl{xe- d b) \ Ms>) ) \ Kb) J / _0 , 

T(l + a + b) \ a/ 



Using (P.10| ) we trade the 1/a factor for the logarithmic integral of the exponent of the pt 
radiator with the rescaled argument, as we did in ( |B.11|) , to obtain 



K ) Ji y r(i + n'(x) + w 



J x z r(l + K' + K'(z)) 1 ; 

The integrand is monotonically decreasing with z. Extracting the pt distribution, we write 

q- J?' e ~ 2n(x) H( ) H(r] _ f°° dz K(x)-TKz) T ( 1 + m ') / D ox 
b ~ ^T(1+2W) {Xh [) ~J X z 6 T{l+K> + K>{z))" [U - 6) 

where we have dropped the factor containing the ratio of A functions as it produces a negli- 
gible effect. Taken together with the regular piece in Q3.24 ) this leads to the expression for 
the shift reported in ( |3.25|) . 

For 1Z' -C 1 we can expand the ratio of the T functions in ( [D . 3| ) to get 

H(x) = / ^ e *«-*M [ 1 - lE R! + lE K'(z) + . . .} 
= (1 - lE Rl) ■ E(x) + lE + O (a s E(x)) . 
The answer for H(B^ 1 ) depends on the relation between 1Z' and y/TZ" ~ •y/cv 
For s < 1 (W < y/cQ ( FXID gives 

= + 7E + o W = (me + ^ + 1 - + 1E 

where we have substituted, according to ( |2.41| ), B ~ Be~ lE (for small 72.' one has \(1Z')/2 = 
1 + 0{K')). 

For s ^> 1 (7V ^> -y/o^) we have instead 

HiB- 1 ) ~ ft' -1 , 

which contribution to the shift becomes negligible, 

if 1 1 



\D T \ W-\hiB + const] asln 2 ^ 

35 



~ s -2 < 1 



For TZ' ^> 1, the contribution of H is even smaller, H(B 1 ) -C (1Z') 1 , due to additional 
fall-off of the T-f unction in the integrand of ( p.3| ). 

We conclude that the approximate expression 0D.4T) , 

H(B- V ) ~ (1 - lE W) EiB- 1 ) + lE , (D.6) 
with B given in ( |2.41| ), can be used everywhere, for arbitrary lni? values. 



E Numerical cross-checks and illustrations 



In this Appendix we present a numerical approach which allows us to illustrate some im- 
portant properties of the results and to estimate (some) next-to-next-to leading effects not 
accounted for by the approximate analytical expressions derived in the paper. 

We express the factorised n-soft-parton cross section through an equation governing a 
branching process, which can be implemented as a Monte Carlo "event" generator — where 
by event one means an ensemble of soft gluons. For the accuracy to which we consider the 
jet broadening, it will be sufficient to consider just the emissions of gluons from the quark 



and antiquark, ignoring their subsequent branching [12 



Considering for the time being only a single jet, if the last gluon was produced with trans- 
verse momentum fcj_i, the probability distribution for the next gluon's transverse momentum 
(h,i < h,i-i) is given by 



k 



dP 



2a s 



ln ]~~ ~ \ ) -M/,7.,_j. /,-,., ' • 



dktj \ iH,i 

where a s = ol s CfI^ and A resums the virtual corrections: 

/■fct.i-i / 
In A(k t ,i-i, h,i) = - -r- 
Jkt ti k t 

Technically, the branching starts from a "fake" gluon, i — 0, with 



(E.l) 



h 4 



k 



f.O 



Qe 



-3/4 



(E.2) 



(E.3) 



which does not contribute to any of the momentum sums. The azimuthal direction of each 
kn is chosen at random. The quark recoil, pt is 



Pt 



8=1 



(E.4) 



and the broadening, as in (|2.1[ ), is 



B 



2^(p* + E|^ 



(E.5) 
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Implementing the branching as a Monte Carlo generator makes the the vector sum extremely 
simple — and correspondingly so the analysis of the mean lnp t as a function of B, or also 
averaged over B. 



If one uses a s (k t ^) in the CMW, or "physical" scheme [ 13 1 , then apart from the overall 



normalisation, the accuracy of the resulting description is the same as that of the resummed 



expressions of [12,24], namely next-to-leading-logarithmic in the exponent. 



In practice, when comparing with the analytical results given in the main part of the 
paper it is most informative to consider the results in the limit of fixed a s (i.e. setting 
/3q = 0) so as to avoid having to introduce an arbitrary infrared cutoff in the numerical 
calculation. 

We recall that the power correction to the broadenings depends linearly on the the 
logarithm of the quark transverse momentum in each jet. Accordingly, from (|3.20| ) we 
expect that (lnp t /Q) — In B in a single jet with broadening B, should be a function only of 
R'(B), and we test this by plotting (\iap t /Q) —\nB versus R'. The resulting curves should 
then be independent of a s . Figure |]a shows such curves for two values of a s , compared with 
the theoretical result. In the Monte Carlo results, there is a small dependence on the value 
of a s — corresponding roughly to a shift of the curves by an amount of order a s . This is a 
sublogarithmic effect and so beyond the accuracy of the analytical calculations (and in any 
case beyond the predictive value of the Monte Carlo). For small values of a s there is good 
agreement between the Monte Carlo and the analytical results. 

When considering the total jet broadening, B T , we just generate two independent single- 
jet configurations (|2.3|) , and construct the combination (lnp^/Q + lnp^/Q)— 2 \\iBt which 
is a function of R'(Bt) only for large values of R'. For small R', it goes as 1/ y/a^. Accordingly, 
in figure ^b the curves for different values of a s differ at smaller R'. However the analytic 
predictions and the Monte Carlo results remain in good agreement, except roughly at the 
level of a shift of order a s as in the case of the single jet curves. 

While the above figures demonstrate the agreement between the Monte Carlo and the 
analytical results, they do not illustrate the shifts themselves as a function of B. For the 
case of a fixed coupling, this is done in figure |5|. The main features are the following: D\y is 
practically independent of a s and almost equals In 1/B. For very small B, is practically 
twice In 1/B. For larger B one can see that the slope of the Dt curves tends to that of the 
D\y curves. Finally at large B one sees an offset in Dt which increases as 1/ Jot^. 



F Collection of final formulas 



We collect here for convenience the final expressions for the broadening distributions and 
means, which include 1/Q confinement effects and were used for the phenomenological anal- 
ysis presented in section [|. 
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Theory, a s = 0.1 
MC, a s = 0.1 
Theory, oc s = 0.01 
MC, oc s = 0.01 



Figure 4: (a) (lnp t /Q) — InBi as a function of 1Z'. The Monte Carlo results are shown for 
various values of a s while the analytic curve is independent of a s ; (b) (lnp^/Q + Inpt^/Q) — 
2 In Bt from the Monte Carlo and from ( |3.25| ) . The Monte Carlo and theory curves for 
a s = 0.01 lie on top of each other. 
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F.l Expressions for the shifts 

The integrated wide-jet broadening distribution: 

-B 



o- 1 I dB w -^- = E W (B) = y!? T \B-\VD 1 {B)) (F.l) 
Jo an w 

with the -B-dependent shift D\ given byQ 

d 1 (b) = inB- 1 + Vo -2-p(n') + x(n') + tp(i + n')-4j(i), 

n , = 2Cf <BQ}_ , lnS _! _ ax -0.6137056 , ' 2) 

7T v 7 

and ip(z) the derivative of the logarithm of T(z). The functions p and x are 



o 



[A(a)]- a = / dz 



o 



1 + z 
2z 



(F.3) 



The integrated total broadening distribution: 

a ~ X I* dB TTET = = V [ t T \B-\VD t {B)) (F.4) 

3 We remark that in 72' the 3/4 is beyond the accuracy that we control, however we choose to keep it 
because it has a clear origin, and is among the next-to-leading corrections to the power-suppressed contri- 
bution. 
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with the -B-dependent shift given by 



D T {B) = 2D 1 (B) + 2[iP(l + 2K')-ip(l + K')] + H(B 



H(x) 



dz 



r(i + 2W) 



t(i + w + n'(z)) ' 



B 



2B 



(F.5) 



where zq corresponds to the position of the Landau pole in the perturbative radiator TZ(z), 
where the integrand vanishes. The form that we use for TZ is the two-loop radiator with the 
one-loop coupling, which, in the physical (cmw) scheme has the simple expression 



TZ(x) 



4Ci 



L - 



In 1 



In a; 



+ Inx 



(F.6) 



where L = InQ/A = 2ir/ (j3 acMw(Q)), and the physical coupling ocmw is related to the 
standard a^g by 



«CMW = «MS [1 + K 



MS 



2tt 



(F.7) 



with 



K = Cz 



67 
18 



7T" 

"6 



-n 



9" f: 3 3 ' 

The two-loop radiator with the two-loop coupling is given in Appendix [A.l| in (|A.23|)P|. 



11C A 2n 



f 



(F.8) 



F.2 PT spectra 



The resummed expressions Q2.42p that we derived for the pt spectra are applicable in the 



small-I? region, and have next-to-leading-logarithmic accuracy. For the purposes of phe- 
nomenology, it is necessary to extend the domain of validity of the perturbative spectra 
towards larger values of B, where it is sufficient to use a fixed calculation. Procedures (the 
R and log-i? matching schemes) for combining the fixed order and resummed results are 
explained in detail in for the thrust and heavy-jet mass distributions, and are directly 
applicable also to the case of the broadenings. The requirement for the use of these proce- 
dures (in particular i?-matching) is that the resummed perturbative answer have only the 
following terms 

oo oo 

In E(V) = da s + Gn, n+ ia r : ln n+1 V + ^ G n , n a n s In" V , (F.9) 

n=l n=l 

and it mustn't have terms such as In V. Suitable expressions for the broadenings were 



presented in ||24j| , equations (18-22). As was shown in |12[ these answers have to be modified 



by an additional factor (both for the wide and total broadenings): 

27?.' 



E(PT,(B) = (jm) ' s - i(B) - (F ' 10) 



'For enthusiasts only! 
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It is vital that 7Z' here be taken as: 



K = 2C " (Q)Cf ^ . (F.11) 

We stress that the 3/4 in 7V is a next-to-next-to-leading effect, and as such is taken care of 
by the matching procedure^. For that procedure to remain intact, it must not be included 
in (TO) . 



F.3 Means 

The leading power correction to the mean total broadening is 

(B) T - (B)™ = V ( 71 + ^-J^ + r]0 + O (yffi) ■ (F-12) 

\2^C F a C Mw(Q) 4 Q °f J 

and the correction to the mean wide-jet broadening is 

<B) »- - ^ - I (w + i " + * + ° ( H ' (F ' 13) 

Here, Q = Qe~ 3 / 4 . The use of Q rather than Q as the scale for ot s , and the choice of «cmw 
rather than ajjg both affect the results at the level of a O (y/c^) term, which formally we 
do not control. However we prefer to keep these corrections since they have clear physical 
origins (the e~ 3//4 factor in the scale has about a 5% effect on the fitted value of ao, while 
the change from cmw to ms schemes has much less effect). 



F.4 Non-PT parameter 

In order to accurately define the non-perturbative parameter V the problem of merging the 
pt and np contributions should be addressed. The relevant procedure was discussed in detail 
in 0. It includes introducing an infrared matching scale \ii (typically chosen to be /if =2 
GeV) and the non-perturbative /independent phenomenological parameter a Q1.3Q which 
quantifies the intensity of QCD interaction over the infrared domain, k < fj,j. 

An explicit expression for V depends on the order to which the perturbative contribution 
is computed, as well as on the scheme. At two-loop level, in the ms scheme, we have 

V = ^-^f {«o(/i/)-« s -/3 ^(ln£ + ^ + l)} ; a a = a m (Q). (F.14) 

The term proportional to K accounts for mismatch between the Ms and the physical scheme, 
with K given above in ( [F.8Q . M. in ( |F.14| ) is the Milan factor resulting from the two- loop 



3 We note that the (3q that we use in this paper differs from that in B4j by a factor of 47r. 
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analysis discussed in Appendix |A.2| (see also |J). This factor is universal for all 1/Q jet 
observables considered in e + e~ annihilation and DIS processes [|23| and reac 



M 



1 + (3a 1 (1.575CU ~ 0.104n/) = 1.490 (1.430) for n f = 3 (0) . 



(F.15) 



The perturbative terms in (|F.14|) (proportional to a s and al) represent the start of the 
series responsible for subtracting off the infrared renormalon divergence in the perturbative 
contribution to the observable. 
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